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1.  Introduction 


This  report  documents  a  set  of  functions,  written  in  C++,  that  can  be  used  to  create,  position,  and 
rotate  rectangles,  as  well  as  test  for  intersection  between  rectangles.  The  functions  build  on 
functions  from  the  y2DOps  namespace.1  The  functions  have  been  grouped  into  the 
yRectangle2D  namespace,  which  is  summarized  at  the  end  of  this  report. 


2.  Rectangles  in  Two  Dimensions 


Suppose  that  a  rectangle  R  is  defined  by  four  position  vectors  R0 ,  Rt ,  R2,  and  R3  as  shown  in 
figure  1. 
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Figure  1 .  A  rectangle  with  comers  defined  by  position  vectors  R0 ,  R{ ,  R2 ,  and  R3 . 
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Laboratory:  Aberdeen  Proving  Ground,  MD,  June  2013. 
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The  functions  described  in  this  report  operate  on  rectangles  that  are  stored  in  eight-element 
arrays,  where  an  array,  R,  is  defined  such  that 


^  ’  Rq ,y  ’^l,x’  ^1  ,y  >^2,x>  ^ 2,y  ’  ^3,x  ’  ^lr  }  • 


(1) 


3.  Creating  Rectangles  -  the  NewRectangle2D()  Function 


The  NewRectangle2D()  function  can  be  used  to  create  rectangles  that  are  placed  with  their 
lower-left  comers  at  the  origin,  as  shown  in  figure  2. 
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Figure  2.  A  rectangle  created  using  the  NewRectangle2D()  function. 

Note  that  the  NewRectangle2D()  function  uses  the  “new”  command  to  allocate  memory  for  the 
array  that  is  pointed  to  by  its  return  value.  Thus,  to  avoid  memory  leaks,  each  use  of  the 
NewRectangle2D()  function  should  be  accompanied  by  a  use  of  the  “delete[]”  operator. 

3.1  NewRectangle2D()  Code 

inline  double*  NewRectangle2D(//<======================CREATES  A  NEW  RECTANGLE 

double  w ,//< - - - - - - - WIDTH  OF  RECTANGLE 

double  l){//< - - - - - - - LENGTH  OF  RECTANGLE 

double  *R=new  double[8];/*<-*/R[0]=R[l]=R[3]=R[6]=0,R[2]=R[4]=w,R[5]=R[7]=l; 
return  R; 

} ;  /  / - YAG  E  NAUT@GMAI  L .  COM  f\j  t\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  LAST~UPDATED~10JUL2013  t\j  r\j  r\j  r\j 


3.2  NewRectangle2D()  Parameters 

w  w  specifies  the  width  of  a  rectangle. 
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1  specifies  the  length  of  a  rectangle. 


3.3  NewRectangle2D()  Return  Value 

The  NewRectangle2D()  function  returns  a  pointer  to  a  new  rectangle. 

3.4  NewRectangle2D()  Example 

The  following  example  begins  by  using  the  NewRectangle2D()  function  to  create  a  rectangle  that 
is  2.0  units  wide  and  1.0  unit  long.  The  coordinates  of  the  comers  of  the  rectangle  are  then 
printed  to  the  screen.  Finally,  the  memory  that  was  allocated  using  the  NewRectangle2D() 
function  is  deallocated  using  the  delete  []  command. 


ftinclude  <cstdio>// . printfQ 

#include  "y_rectangle_2d . h" 
int  main(){ 

double*  R=yRectangle2D: :NewRectangle2D(2jl); 
printf("  x  ,  y\n"); 

for(int  i=0; i<4;++i)printf ("  %4.1f  ,  %4.1f\n"JR[2*i],R[2*i+l]); 
delete[]  R; 

>;//  f\J  f\J  f\J  f\J  f\J  V  AG  E  NAUTfgJGMAI  L .  LAST  ~U  P  DAT  E  D~  1 0  D  U  L2013~~~~~~ 


OUTPUT: 


X 

y 

y 

0.0 

y 

0.0 

2.0 

y 

0.0 

2.0 

y 

1.0 

0.0 

y 

1.0 

4.  Translating  Rectangles  -  the  TranslateRectangle2D()  Function 


The  TranslateRectangle2D()  function  can  be  used  to  reposition  rectangles  (without  rotation)  by  a 
displacement  vector  d ,  as  shown  in  figure  3. 
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Figure  3.  A  rectangle  that  has  been  repositioned  using  the  TranslateRectangle2D()  function. 

4.1  TranslateRectangle2D()  Code 

inline  void  TnanslateRectangle2D(//<====================TRANSLATES  A  RECTANGLE 

double  R[8]j//< - - - A  RECTANGLE  (CREATE  USING  NewRectangle2D( ) ) 

const  double  d[2]  ){//< - - - - - DISPLACEMENT  VECTOR 

fon(int  i=0; i<4;++i)y2D0ps : :Translate2D(R+2*i,d); 

};// - YAGENAUT@GMAIL.COM - LAST~UPDATED~10DUL2013 - 


4.2  TranslateRectangle2D()  Parameters 

R  R  specifies  a  rectangle  that  has  been  created  using  the  NewRectangle2D()  function. 

d  d  specifies  a  two-element  array  that  stores  the  displacement  vector  d  (d=  {dx ,  dy } ).  d 

determines  the  magnitude  and  direction  by  which  R  is  translated. 

4.3  TranslateRectangle2D()  Example 

The  following  example  begins  by  using  the  NewRectangle2D()  function  to  create  a  rectangle  that 
is  2.0  units  wide  and  1.0  unit  long.  The  rectangle  is  then  translated  by  3.0  units  in  the  x 
direction  and  2.5  units  in  the  y  direction  using  the  TranslateRectangle2D()  function.  The 

coordinates  of  the  comers  of  the  translated  rectangle  are  then  printed  to  the  screen.  Finally,  the 
memory  that  was  allocated  using  the  NewRectangle2D()  function  is  deallocated  using  the 
delete[]  command. 
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#include  <cstdio>// . printf() 

#include  "y_rectangle_2d . h" 
int  main(){ 

double*  R=yRectangle2D: :NewRectangle2D(2jl); 
double  d[2]={3j2.5}; 

yRectangle2D: :TranslateRectangle2D(R.,d) ; 
printf("  x  y\n"); 

for(int  i=0; i<4;++i)printf ( "  %4.1f  ,  %4.1f\n",R[2*i],R[2*i+l]); 
delete[]  R; 

};//  YAGENAUT@GMAIL.COM  f\J  f\J  LAST~UPDATED~101UL2013 


OUTPUT: 


X 

y 

y 

3.0 

y 

2.5 

5.0 

y 

2.5 

5.0 

y 

3.5 

3.0 

y 

3.5 

5.  Rotating  Rectangles  -  the  RotateRectangle2D()  Function 


The  RotateRectangle2D()  function  can  be  used  to  rotate  rectangles  about  their  centers  by  an 
angle,  0 ,  as  shown  in  figure  4. 


y 


Figure  4.  A  rectangle  that  has  been  rotated  using  the  RotateRectangle2D()  function. 


5 


5.1  RotateRectangle2D()  Code 


inline  void  RotateRectangle2D(//<=========R0TATES  A  RECTANGLE  ABOUT  ITS  CENTER 

double  R[8]j//< - - - A  RECTANGLE  (CREATE  USING  NewRectangle2D( ) ) 

double  rads){//< - - - THE  ANGLE  OF  THE  ROTATION  (CCW  IS  POSITIVE) 

double  M[4]Jo[2]={(R[0]+R[2]+R[4]+R[6])/4,(R[l]+R[3]+R[5]+R[7])/4}; 
y2DOps: :RMatrix2D(M,rads); 

for(int  j=0; j<4;++j)y2D0ps: :Rotate2D(R+2*j,o,M); 

};// - YAGENAUT@GMAIL.COM - LAST~UPDATED~103UL2@13 - 


5.2  RotateRectangle2D()  Parameters 

R  R  specifies  a  rectangle  that  has  been  created  using  the  NewRectangle2D()  function. 

rads  rads  specifies  the  angle  (in  radians)  of  a  rotation.  The  direction  of  the  rotation  is 
counterclockwise. 

5.3  RotateRectangle2D()  Example 

The  following  example  begins  by  using  the  NewRectangle2D()  function  to  create  a  rectangle  that 
is  2.0  units  wide  and  1.0  unit  long.  The  rectangle  is  then  rotated  by  it  1 2  using  the 
RotateRectangle2D()  function.  The  coordinates  of  the  comers  of  the  rotated  rectangle  are  then 
printed  to  the  screen.  Finally,  the  memory  that  was  allocated  using  the  NewRectangle2D() 
function  is  deallocated  using  the  delete  []  command. 

ftinclude  <cstdio>// . printfQ 

ftinclude  "y_rectangle_2d . h" 
int  main(){ 

double  w=2,l=l;// . size  of  rectangle 

double*  R=yRectangle2D: :NewRectangle2D(w,l); 
yRectangle2D: : RotateRectangle2D(R, 3. 14159265358979/2); 
printff'  x  ,  y\n")j 

for(int  i=0; i<4;++i)printf ( "  %4.1f  ,  %4.1f\n"JR[2*i],R[2*i+l]); 
delete[]  R; 

};//  r\j  r\j  r\j  r\j  r\j  YAGENAUT@GMAIL.COM  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  LAST~UPDATED~10JUL2013  f\j  f\j  f\j  f\j  f\j 


OUTPUT: 


x  ,  y 
1.5  ,  -0.5 
1.5  ,  1.5 

0.5  ,  1.5 

0.5  ,  -0.5 
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6.  Checking  for  Rectangle  Overlap  -  the  RectangleIntersect2D()  Function 


The  RectangleIntersect2D()  function  can  be  used  to  determine  whether  or  not  any  of  the  sides  of 
two  rectangles  intersect.  Figure  5  shows  examples  of  intersecting  and  nonintersecting  rectangles. 


y 


Figure  5.  Intersecting  rectangles:  A  and  B  intersect,  A  and  C  do  not. 


6.1  RectangleIntersect2D()  Code 

inline  bool  RectangleIntersect2D(//<===D0  THE  SIDES  OF  2  RECTANGLES  INTERSECT? 


const  double  Rl[8],//< - RECTANGLE  #1  (CREATE  USING  NewRectangle2D( ) ) 

const  double  R2[8],//< - RECTANGLE  #2  (CREATE  USING  NewRectangle2D( ) ) 


double  e=lE-9){//<- -CUTOFF  VALUE  FOR  DETERMINING  IF  TWO  SIDES  ARE  PARALLEL 
for(int  j=0; j<4;++j)for(int  j2=0; j2<4;++j2){ 
double  t[2],x[2]; 

double  LA[4]={Rl[2*j  ],Rl[2*j  +1]  ,  Rl[(j  +1)%4*2] , Rl[ ( j  +1)%4*2+1]}; 
double  LB[4]={R2[2*j2] , R2[2*j2+1]  ,  R2[ ( j2+l)%4*2] , R2[ ( j2+l)%4*2+l] }; 
if ( !y2D0ps : : IParameters2D(t, LA, LB) ) continue; 
if(y2D0ps: : Intersect2D(x,t, LA) )return  true;} 
return  false; 

} ;  /  / - YAG  E  NAUT@GMAI  L .  COM  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  LAST~UPDATED~10DUL2013  f\j  f\j  f\j  f\j  f\j  f\j 


6.2  RectangleIntersect2D()  Parameters 

R1  R1  specifies  a  rectangle  that  has  been  created  using  the  NewRectangle2D()  function. 

R2  R2  specifies  a  second  rectangle  that  has  been  created  using  the  NewRectangle2D() 
function. 

e  e  specifies  the  cutoff  value  for  determining  if  two  sides  are  parallel. 


7 


6.3  RectangleIntersect2D()  Return  Value 

The  RectangleIntersect2D()  function  returns  true  if  any  side  of  R1  intersects  any  side  of  R2, 
except  for  the  special  case  where  two  sides  intersect  at  more  than  one  point. 

6.4  RectangleIntersect2D()  Example 

The  following  example  begins  by  creating  and  positioning  the  three  rectangles  shown  in  figure  5. 
The  RectangleIntersect2D()  function  is  then  used  to  determine  whether  or  not  the  rectangles 
intersect.  The  results  are  printed  to  the  screen.  Finally,  the  memory  that  was  allocated  using  the 
NewRectangle2D()  function  is  deallocated  using  the  delete[]  command. 


#include  <cstdio>// . printfQ 

#include  "y_rectangle_2d . h" 
int  main(){ 

double*  A=yRectangle2D: :NewRectangle2D(2jl); 
double*  B=yRectangle2D: :NewRectangle2D(2.,l); 
double*  C=yRectangle2D: :NewRectangle2D(6.,3); 
double  d[2]={lj . 5); 

yRectangle2D: :TranslateRectangle2D(B,d); 
d[0]=- . 5,d[l]=- . 5; 

yRectangle2D: :TranslateRectangle2D(C,d); 
printf ("A  intersects  B?  %s\n", 

yRectangle2D: : RectangleIntersect2D(Aj  B) ?"true" : "false"); 
printf("B  intersects  C?  %s\n", 

yRectangle2D: : RectangleIntersect2D(BjC) ?"true" : "false"); 
printf ("C  intersects  A?  %s\n", 

yRectangle2D: : RectangleIntersect2D(C.,  A)  ?"true" : "false"); 
delete[]  A; 
delete[]  B; 
delete[]  C; 

};//  YAGENAUT@GMAIL.COM  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  r\j  LAST~UPDATED~10JUL2013  f\j  f\j  f\j  f\j  f\j  f\j 


OUTPUT: 


A  intersects  B?  true 
B  intersects  C?  false 
C  intersects  A?  false 
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7.  Interior/Exterior  Check  -  the  RectangleInterior2D()  Function 


The  RectangleInterior2D()  function  can  be  used  to  determine  whether  or  not  a  point  lies  inside  a 
rectangle. 
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Figure  6.  Interior  Points:  P  is  interior  to  B  ,  but  not  to  A  . 

7.1  RectangleInterior2D()  Code 


inline  bool  RectangleInterior2D(//<=============IS  A  POINT  INSIDE  A  RECTANGLE? 

const  double  R[8],//< - A  RECTANGLE  (CREATE  USING  NewRectangle2D( ) ) 

const  double  P[2]){//< - - - - - - - A  POINT 

double  A[4]={P[0],P[1]jP[0]+R[6]-R[0]jP[1]+R[7]-R[1]}; 
double  B[4]={P[0],P[1]JP[0]+R[2]-R[0],P[1]+R[3]-R[1]}; 
double  C[4]={R[2],R[3]jR[4],R[5]}; 
double  D[4]={R[6],R[7]jR[4],R[5]}; 
double  tl[2];/*<-*/y2D0ps: :IParameters2D(tl,AJD); 
double  t2[2];/*<-*/y2D0ps: :IParameters2D(t2,BJC); 
return  tl[0]>0&&tl[0]<l&&t2[0]>0&&t2[0]<l; 

};// - YAGENAUT@GMAIL.COM - LAST~UPDATED~10DUL2013 - 


7.2  RectangleInterior2D()  Parameters 

R  R  specifies  a  rectangle  that  has  been  created  using  the  NewRectangle2D()  function. 

P  P  specifies  a  point  that  lies  in  the  plane  of  R. 

7.2  RectangleInterior2D()  Return  Value 

The  RectangleInterior2D()  function  returns  true  if  P  lies  interior  to  R,  and  false  otherwise. 
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7.2  RectangleInterior2D()  Example 

The  following  example  begins  by  creating  and  positioning  the  rectangles  and  the  point  that  are 
shown  in  figure  6.  The  RectangleInterior2D()  function  is  then  used  to  determine  whether  or  not 
the  point  is  interior  to  either  of  the  rectangles.  The  results  are  printed  to  the  screen.  Finally,  the 
memory  that  was  allocated  using  the  NewRectangle2D()  function  is  deallocated  using  the 
delete[]  command. 

ftinclude  <cstdio>// . printfQ 

ftinclude  "y_rectangle_2d . h" 
int  main(){ 

double*  A=yRectangle2D: :NewRectangle2D(2,l); 
double*  B=yRectangle2D: :NewRectangle2D(2jl); 
double  d[2]={3j2); 

yRectangle2D: :TranslateRectangle2D(B,d); 

double  P[2]={4,2.5}; 

printf ("P  interior  to  A?  %s\n", 

yRectangle2D: : RectangleInterior2D(A, P) ?"true" : "false") ; 
printf("P  interior  to  B?  %s\n", 

yRectangle2D: : RectangleInterior2D(Bj  P) ?"true" : "false"); 
delete[]  A; 
delete[]  B; 

>5// -  — YAG  E  NAUTfaJGMAI  L .  LAST  ~U  P  DAT  ED~10DUL2013~~~~~~ 


OUTPUT: 


P  interior  to  A?  false 
P  interior  to  B?  true 


8.  Example  -  Placing  TV  Uniformly  Distributed  Rectangles  in  a  Rectangular 
Region 


The  following  example  uses  functions  from  the  yRectangle2D  namespace  to  place  800  rectangles 
that  are  2.0  units  by  3.0  units  in  size  within  a  200.0  unit  by  100.0  unit  rectangle.  Figure  6 
provides  a  visual  of  the  800  rectangles.  Note  that  none  of  the  rectangles  intersect,  either  with 
each  other,  or  with  the  bordering  rectangle.  The  example  then  creates  1,000,000  randomly 
positioned  points  and  tests  to  see  if  any  of  the  points  are  interior  to  any  of  the  rectangles. 

Figure  7  provides  a  visual  of  the  points  that  are  interior  to  rectangles. 
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#include  <cstdio>// . FILE f reopen ( ) , stdout.,  f  close  ( ) , printf ( ) 

#include  <cstdlib>// . rand ( ) , RAND_MAX 

#include  "y_rectangle_2d . h" 
int  main(){ 

int  N=800;// . number  of  rectangles  in  region 

double  w=2.,l=3;// . size  of  rectangles 

double  1aI=200,  L=100;  / / . size  of  region 

double*  B=yRectangle2D:  :NewRectangle2D(WJ  L);// . region  border 

double**  R=new  double*[N];// . storage  for  a  set  of  rectangles 

FILE  *f=freopen(Mrectangles.txt,V,w",stdout);// . redirect  output  to  a  file 

for(int  i=0; i<N;++i){ 

R[i]=yRectangle2D: :NewRectangle2D(w, 1); 

yRectangle2D: :RotateRectangle2D(R[i],2*3.14159265358979*rand()/RAND_MAX); 

redo:// . if  two  rectangles  intersect,  come  back  to  this  point 

double  d[2]={-w/2+W*rand()/RAND_MAX,  -l/2+L*rand( )/RAND_MAX}; 
yRectangle2D: :TranslateRectangle2D(R[i] ,d); 
if (yRectangle2D: : RectangleIntersect2D(R[i] , B) ){ 

d[0]*=-l,d[l]*=-l,yRectangle2D: :TranslateRectangle2D(R[i] ,d); 
goto  redo;} 

for (int  i2=0; i2<i;++i2)if (yRectangle2D: : RectangleIntersect2D(R[i] , R[i2] ) ){ 
d[0]*=-l,d[l]*=-l,yRectangle2D: :TranslateRectangle2D(R[i] ,d); 
goto  redo;} 

printf(''%f,%f,%f,%f,%f,%f,%f,%f\n%R[i][0],R[i][l],R[i][2],R[i][3], 

R[i][4],R[i][5],R[i][6],R[i][7]);} 


fclose(f); 

freopen("points.txt"J  "w'^stdout);// . redirect  output  to  a  file 

for(int  k=0;k<1000000;++k){ 

double  P[2]={W*rand()/RAND_MAX,L*rand()/RAND_MAX}; 
for (int  i=0; i<N;++i)if (yRectangle2D: : RectangleInterior2D(R[i] , P) ){ 
printf ("%fJ%f\n",P[0]JP[l]); 
break; }} 

for(int  i=0;i<N;++i)delete[]  R[i]; 
fclose(f); 

};//  YAGENAUT@GMAIL.COM  C\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  f\J  DRAFT  t\j  f\j  f\j  f\j  f\j  f\j  f\j  f\j  rsj  LAST~UPDATED~10DUL2013  f\j  t\j  t\j  f\j  t\j  t\j 
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Figure  7.  N uniformly  distributed  rectangles. 


Figure  8.  Randomly  generated  points  that  are  interior  to  800  uniformly  distributed  rectangles. 
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9.  Summary 


A  summary  sheet  is  provided  at  the  end  of  this  report.  It  presents  the  yRectangle2D  namespace, 
which  contains  the  five  functions  that  are  described  in  this  report.  Also  presented  is  the  y2DOps 
namespace,  which  is  utilized  by  functions  in  the  yRectangle2D  namespace,  and  the  example  that 
is  presented  in  section  8. 
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y_  rectangle_2d.h 


#ifndef  Y_RECTANGLE_2D_GUARD 
#define  Y_RECTANGLE_2D_GUARD 
#include  "y_2d_ops.h" 
namespace  yRectangle2D{ 

inline  double*  N<j»wRectaMM(7/<======================CREATES  A  NEW  RECTANGLE 

double  w ,//< - WIDTH  OF  RECTANGLE 

double  l){//< - LENGTH  OF  RECTANGLE 

double  *R=new  double[8] j/*<-*/R[0]=R[l]=R[3]=R[6]=0JR[2]=R[4]=wJR[5]=R[7]=l; 
return  R; 

} ;  // - YAGENAUT@GMAI  L .  COM - - LAST~UPDATED~10DUL2013 - 

inline  void  TranslateRectangle2D(//<====================TRANSLATES  A  RECTANGLE 

double  R[8] ,/ / < - A  RECTANGLE  (CREATE  USING  NewRectangle2D() ) 

const  double  d[2]){//< - DISPLACEMENT  VECTOR 

for(int  i=0;i<4;++i)y2DOps: :Translate2D(R+2*ijd); 

} ;  // - YAGENAUT@GMAI  L .  COM - - LAST~UPDATED~10DUL2013 - 

inline  void  ( / /<=========rotates  a  rectangle  about  its  center 

double  R[8] ,/ /< - A  RECTANGLE  (CREATE  USING  NewRectangle2D()  ) 

double  rads){//< - THE  ANGLE  OF  THE  ROTATION  (CCW  IS  POSITIVE) 

double  M[4]jO[2]={(R[0]+R[2]+R[4]+R[6])/4j(R[1]+R[3]+R[5]+R[7])/4}; 
y2DOps : : RMatrix2D(M,  rads) ; 

for(int  j=0; j<4;++j)y2D0ps: : Rotate2D(R+2*jJoJM); 

} ;  // - YAGENAUT@GMAI  L .  COM - - LAST~UPDATED~10DUL2013 - 

inline  bool  RectangleIntersect2D(//<===DO  THE  SIDES  OF  2  RECTANGLES  INTERSECT? 

const  double  Rl[8],//< - RECTANGLE  #1  (CREATE  USING  NewRectangle2D() ) 

const  double  R2[8]J//< - RECTANGLE  #2  (CREATE  USING  NewRectangle2D() ) 

double  e=lE-9){/ / <- -CUTOFF  VALUE  FOR  DETERMINING  IF  TWO  SIDES  ARE  PARALLEL 
for(int  j=0; j<4;++j)for(int  j 2=0; j2<4;++j2){ 
double  t[2]jX[2]; 

double  LA[4]={Rl[2*j  ]JRl[2*j  +1]  ,  Rl[(j  +1)%4*2] , Rl[ ( j  +1)%4*2+1]>; 
double  LB[4]  =  {R2[2*j2]j  R2[2*j2+1]  ,  R2[ ( j2+l)%4*2] , R2[ ( j2+l)%4*2+l]}; 
if ( !y2DOps: : IParameters2D(t ,  LA, LB) )continue; 
if(y2DOps: : Intersect2D(xJtJ LA))return  true;} 
return  false; 

};// - YAGENAUT@GMAIL.COM - LAST~UPDATED~10DUL2013 - 

inline  bool  RectangleInterior2|(//<=============IS  A  POINT  INSIDE  A  RECTANGLE? 

const  double  R[8]J//< - A  RECTANGLE  (CREATE  USING  NewRectangle2D() ) 

const  double  P[2]){//< - A  POINT 

double  A[4]={P[0]jP[1]jP[0]+R[6]-R[0]jP[1]+R[7]-R[1]}; 

double  B[4]={P[0]jP[1]jP[0]+R[2]-R[0]jP[1]+R[3]-R[1]}; 

double  C[4]={R[2]jR[3]jR[4]jR[5]}; 

double  D[4]={R[6]jR[7]jR[4]jR[5]}; 

double  tl[2] ;/*<-*/y2D0ps : : IParameters2D(tljAjD) ; 

double  t2[2] ;/*<-*/y2D0ps : : IParameters2D(t2j  BjC) ; 

return  tl[0]>0&&tl[0]<l&&t2[0]>0&&t2[0]<l; 

} ;  // - YAGENAUT@GMAIL .  COM - LAST~UPDATED~10DUL2013 - 

}// @@@@@@@@@@@@@@@@@@@@@@@@@@@( 

#endif 
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EXAMPLE 


#define  main  yRectangle2DRectangularBattlefieldExample//@@@@@@@@@@@@@@@@@@@@@@@@ 

#in elude  <cstdio>// . FILE, freopen( ) , stdout, f close () jprintf() 

#include  <cstdlib>// . rand( )  ,  RAND_MAX 

#include  "y_rectangle_2d. h" 
int  main(){ 

int  N=800;// . number  of  rectangles  in  region 

double  w=2jl=3;// . size  of  rectangles 

double  W=200j  L=100; / / . size  of  region 

double*  B=yRectangle2D:  :NewRectangle2D(Wj  L);// . region  border 

double**  R=new  double*[N] ;// . storage  for  a  set  of  rectangles 

FILE  *f=f reopen("rectangles .txt ", "w" , stdout) ;// . redirect  output  to  a  file 

for(int  i=0;i<N;++i){ 

R[i]=yRectangle2D : : NewRectangle2D(Wj 1) ; 

yRectangle2D: :RotateRectangle2D(R[i]J2*3.14159265358979*rand()/RAND_MAX); 

redo:// . if  two  rectangles  intersect,  come  back  to  this  point 

double  d[2]={-w/2+W*rand( )/RAND_MAXj -l/2+L*rand( )/RAND_MAX}; 
yRectangle2D: :TranslateRectangle2D(R[i] ,d) ; 
if (yRectangle2D: : RectangleIntersect2D(R[i] ,  B) ){ 

d [0]*=-l, d[l]*=-l,yRectangle2D: :TranslateRectangle2D(R[i] ,d) ; 
goto  redo;} 

for (int  i2=0; i2<i;++i2)if (yRectangle2D: : RectangleIntersect2D(R[i] ,  R[i2] ) ){ 
d[0]*=-l, d[l]*=-l,yRectangle2D: :TranslateRectangle2D(R[i] ,d) ; 
goto  redo;} 

printf ("%f j%f j%f j%f j%f j%f j%f j%f\n" jR[i] [0] JR[i][l]JR[i][2]JR[i][3]J 

R[i][4]JR[i][5]JR[i][6]JR[i][7]);} 


fclose(f) ; 

freopen("points.txt"j  "w"j  stdout);// . redirect  output  to  a  file 

for(int  k=0;k<1000000;++k){ 

double  P[2]={W*rand()/RAND_MAXj  L*rand( )/RAND_MAX}; 
for (int  i=0;i<N;++i)if (yRectangle2D: : RectangleInterior2D(R[i] ,  P) ){ 
printf ("%f ,%f \n",P[0] , P[l] ) ; 
break; }} 

for(int  i=0;i<N;++i)delete[]  R [ i] ; 
fclose(f) ; 

};// - YAGENAUT@GMAIL.COM - DRAFT - LAST-UPDATED-103UL2013 - 
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